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Abstract 

We consider an elliptic billiard whose shape slowly changes. During slow evolution of the billiard 
certain resonance conditions can be fulfilled. We study the phenomena of capture into a resonance 
and scattering on a resonance which lead to the destruction of the adiabatic invariance in the 
system. 
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Billiard systems are important models in different fields of physics and the theory of 
dynamical systems l], Q- El U, S] . Recently the billiards with varying parameters became the 
object of interest |6|. In the present paper we consider effects of scattering on resonances 
and capture into a resonance in the dynamics of a particle in an elliptic billiard with slowly 
varying parameters. We use a scheme of analysis of resonant phenomena in Hamiltonian 
systems possessing slow and fast variables this scheme is a Hamiltonian version of a 
general scheme 0, 0]. These methods were designed and corresponding theorems were 
oroved for smooth systems. In the recent paper on slowly perturbed rectangular billiards 
lfl| it was shown that the methods can be also applied adequately for billiard systems 
which possess discontinuities. The present paper gives another example of a billiard system 
possessing effects of capture into a resonance and scattering on a resonance. 

Consider a particle moving in an elliptic billiard. Let ri,T2 be the particle distances 
from foci (points 0\,Oi correspondingly), and 2c is the distance between the foci. The 
Hamiltonian of the system in elliptic coordinates £ = r\ +r2, tj = r% — r-i has the following 
form [l| : 



p 2 - Ac 2 Ac 2 - n 2 

Ho = W-^f + ^l 2 ^ + m 0) ' (1) 
where {/(£, a) is the potential of the billiard wall: 



. oo, if P > 2a 
U(t,a) = { (2) 
0, if £ < 2a 

We investigate the problem of an elliptic billiard with slowly changing parameters c = 
c(et), a = a(et) which slowly rotates with angular velocity u. Let us start with cartesian 
coordinates (x,y) such that the ellipse foci are located in (c, 0) and (— c, 0). Let (p x ,Py) 
be momenta conjugated to (x,y). Changing to new variables (P u ,u, P v ,v) by means of a 
canonical transformation with a generating function W = c(p x cosh v cos u — p y sinh v sin u) 
(see Ref. we get the Hamiltonian 

P 2 + P 2 ec' P u sin 2u - P v sinh 2v P u sinh 2v + P v sin 2u 

H = -TT, : r H • : V UJ ■ h Li(v, a I C) 

c^(cosh2f — cos 2m) c cosh 2v — cos 2u cosh 2v — cos 2u 



where 



. oo, if cosht> > a/c 
U(v,a/c) = { (4) 
0, if cosht> < a/c 



New coordinates relate to the elliptic coordinates by means of the following formulas: 



Let us introduce the variable P t which is canonically conjugated to time t, and the variable 
t = et. Now we can consider the system as the autonomous three-dimensional system with 
Hamiltonian: H'(P u ,u, P v ,v, P t ,r) = H + P t . Canonically conjugated pairs are (P u ,u), 
(P v , v), and (Pt, This Hamiltonian system is conservative. Let us consider its dynam- 

ics on some level of the Hamiltonian H' = E. Denote k = 2(E — P t ) (in the unperturbed 
system k > 0, since the Hamiltonian (jSJ) is greater than 0). By means of simple canonical 
transformation we can change to Hamiltonian H = H' — E and then introduce new time 
variable t' such that dt/ dt'=c 2 (cosh 2v — cos 2m) 11]. Since the value of the Hamiltonian 
H = H' — E is equal to 0, after introducing new time variable t' new Hamiltonian 7i is 
obtained by multiplying the Hamiltonian H by c 2 (cosh2t> — cos 2m). 
So we get the following Hamiltonian: 

H = Pi - c 2 k cosh 2 v + Pi + c 2 k cos 2 u + U(v, a/c) + ec'c [P u sin 2u - P v sinh 2v] + 



The Hamiltonian © describes the dynamics of a particle in the elliptic billiard perturbed 
by slow deformation and rotation. The unperturbed problem (r=const, e = uj = 0) is 
integrable. In addition to conserved energy there is another constant of motion c\ = c u = c v : 



We introduced variables c u , c v which are equal to each other in the unperturbed system 
and slightly differ in the perturbed system. The unperturbed system can be regarded as 
two uncoupled oscillators whose phase portraits are presented in Fig. 1. One can introduce 
action-angle variables (I u , <f> w I v , 4> v ) by means of canonical transformation with a generating 
function S = S(I U , I v ,u,v, Pt,r) which contains r,P t as parameters. In the new variables 
the unperturbed Hamiltonian JF transforms to 7i = 'Ho(Iu, Iv, Pu t) — 7~Lu(Iu, P, r ) + 
TC V (I V , P*, t). The function S has the form 



£ = 2c cosh v, r] = 2ccosu. 



(5) 



+ cue 2 [P u sinh 2v + P v sin 2u] = JF + eT\ = 0. 



(6) 



P 2 — c 2 k cosh 2 v = c, 
P.. + c k cos u = —c, 



(7) 
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Actions of the system have the following form (see also |12L Il3|): 



Case 1 : c 2 k + c u > 0. 



L = - 



] /-cosh 1 (a/c) 



7T JO 



~\~ kc 



I rcos 

lu = / 

7T Jtt/2 



kc 2 

i _ / — in 



dv\ c v + c 2 k cosh f 



_F(a, m) — V kc 2 E(a, m) + fcau 



a 2 — c 2 



du\/ —c„ — c 2 k cos 2 u = — 



c v + /ca 2 
\Jkc 2 



K(m) + v / fe 2 E(m) 



where m = — 
c 2 k 



a = arcsm 



lk(a 2 -c 2 ) 
fca 2 + c u 



(9) 



Case 2 



"] /-cosh 1 (a/c) 



L, — — I i dv\/c v + c 2 k cosh 2 v 



7T 7 cosh" 1 J 



I u = — [ duJ —c u — c 2 k cos 2 u 
where (3 = arcsin 



\\ k{a 2 - c 2 ) ' 
Frequencies of the system have the form 



dHo 1 2tt v / ^ 



c 2 k + c u < 0. 
1 



7T 



/a 2 /c + c„ 
a 2 — c 2 



^E(i), 
7r m 



dl v F(a,m) 



a;,, 



9c„ 



case 1 
dfto _ J_ 

QT ~ din 
Ulu dc u 



m 



27rV^ 
K(m) ' 



(10) 



UJ r . 



2tt^-c v 



case 2 



2tt v /^c^ 



Now consider the perturbed problem (e ^ 0), following Ref. Q. Let us make in the 
system with Hamiltonian (jUJ) the canonical transformation of the variables 



(P u , u, P v , v, P t , t) -> (I u , (j) u , I v , (j) v , P t , t) 



(11) 



determined by the generating function 



S 2 = P t e t + S(I u ,I v ,u,v,P u t). 



(12) 
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Formulas for the transformation of the variables have the form 

dS 



-, a = u, v, 



P-™ P-™ (13) 

„ = dS dS 
P t = P t + e—, f = r + e- w . 

or oP t 

Hamiltonian © in the new variables is as follows: 

1 Pi (Pui (pvi 

P t ,f) + 0(e 2 ), 
dFodS &H dS 

The variables (I u , (fi u , I v , 4> v , P t , f) are 0(e)- close to the variables (I u , <fi u , I v , (fi v , P t , r). Hence- 
forth the bar over the new variables is omitted and the new Hamiltonian is 

U = H (I U , Pi Pi r) + eH^Iu, Pi <j> U i 4>vi Pi r) + 0(e 2 ). (15) 

The differential equations of the motion have the form 

OT~C cfH. 
I a = -e— — + 0(e 2 ), cj) a = u> a (I u ,I v ,P t ,r) + e— - + 0(e 2 ), a = u,v, (16) 
d(p a dl a 

= ~ £ V + 0(6 >■ T = £ m + £ m + 0(£ >• 

Averaging of the right hand sides of |TT)|) over (p a and discarding terms 0(e 2 ) gives an 
averaged system 

P = 0, P t = -e-*—, t = e— — . (17) 
dr dP t 

Approximation (|T7j) is called an adiabatic approximation [9|. Trajectories of the system (fTTj) 
are called adiabatic trajectories. In the adiabatic approximation P jV = const. The adiabatic 
approximation breaks down in a vicinity of a resonant surface in the phase space where the 
resonance condition k u u u + k v u v = is fulfilled (k u , k v are integers, k 2 + k 2 ^ 0). Although 
resonant surfaces are dense in the phase space of our three-dimensional Hamiltonian system, 
for small e only finite number of low-order resonances are important (the order of a resonance 
is k = \k u \ + \k v \) 0]. 

In the exact (non-averaged) system variables l a are approximate adiabatic invariants, i.e. 
they are well-conserved in a large area of phase space (far from resonant surfaces). Near 
a resonant surface of a given resonance the system (1151) can be transformed into standard 

Q 

"perturbed pendulum-like system" form. Details can be found in Ref. 



During the motion near (k u , k v ) - resonance the magnitude of J = —k v I u + k u I v is 
approximately conserved. Far from resonances the magnitudes of I u , I v are approximately 
conserved. 

So, a phase point in the (I u , I v , P t , r) - space moves in a following way: while it is far from 
low-order resonance surfaces k u u u (I u , P t , r) + k v u v (I v , P t , r) = 0, it moves in a vicinity of an 
adiabatic curve I U)V =const. When it approaches a resonant surface (enters a resonant zone), 
it leaves the adiabatic curve and can be either captured into the resonance or scattered on the 
resonance. In the former case it continues its motion in the vicinity of the resonant surface 
until become ejected from the resonance, whereas in the latter case its adiabatic invariants 
I U} I v undergo jumps ~ y/e and after crossing the resonant zone the particle continues its 
motion along another adiabatic curve. Similar dynamics was investigated recently in Refs. 

li) lh) 3- Corresponding dynamics of the particle in the coordinate space (x, y) can be 
described as follows. In the adiabatic approximation each segment of the particle trajectory 
in a coordinate system rotating with the billiard tangents a confocal quadrics (caustics) 
determined by conditions I u<v = const. The caustics slowly evolves in accordance with 
evolution of the billiard parameters. In the exact system, while a phase point is far from 
low-order resonance surfaces, each segment of the particle trajectory tangents a confocal 
quadrics which is close to the caustics determined by l u<v = const. While captured in a 
(m, n) resonance, each segment of the particle trajectory tangents a confocal quadrics which 
is close to the quadrics determined by conditions J = —nl u + ml v = const, u) u /uj v = —n/m 
= const. 

Consider the case of the billiard whose parameters a, c evolve periodically in time. If 
an adiabatic trajectory crosses a resonant surface, it crosses the surface at the same point 
periodically in time. A phase point which moves near this adiabatic trajectory also crosses 
the resonant surface repeatedly. Accumulation of changes of I u , I v due to multiple passages 
through resonances leads to the destruction of the adiabatic invariance in the system Q. 

It should be noted that the billiard could be deformed in such a way that there will be 
no passages through resonances and the dynamics will be determined by KAM theory 
Let the parameters a and c be changed in such a way that a/c =const. In the system © 
one can change from variables (P t , r) to (P t = c 2 k, f = f(r) = J^ Q ^J) using simple canonical 
transformation which does not change other variables. Then, Hamiltonian Tio depends only 
on actions I U ,I V , Pt and KAM theory could be applied. This means that most of the phase 
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space (/„, I v , P t ) is filled up by invariant tori which are close to the tori I u , I v , P t =const. 

We performed numerical investigations of the system. The results are shown in figures 
2-6. In figures 2-3 the results of numerical investigations of slowly deforming elliptic 
billiards are shown. Semi-axis of a billiard, d a and db, were being changed periodically 
with time. Jumps of the adiabatic invariant I v on a (2,1) resonance are shown in Fig. 2. 
A single jump of the adiabatic invariant is shown in Fig. 3. In Figures 4-6 the results 
of numerical investigations of slowly rotating and deforming billiards are presented. In 
Fig. 4 the capture of a phase point into a (2,1) resonance is shown. The captured point 
moves along a resonant curve until it escapes from the resonance. In Fig. 5 a phase 
point that were initially captured into resonance remains captured forever. Dynamics 
of adiabatic invariant J of that system is shown in Fig. 6. In Fig. 7 we depicted 
results of numerical investigation of a deforming billiard with u> = 0, a/c = const. 
Coordinates (x,y) of a particle at the moments of collisions with a wall are shown. It 
is easy to see that dynamics is regular (compare that picture with Fig. 8, where a/c^ const). 

This work was supported in part by grants RFBR 00-01-00538, RFBR 00-15-96146, and 
INTAS 00-221. A.I. Neishtadt acknowledges support from the program " Integration" -B0053. 
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FIG. 1: Phase portraits of the oscillators (7). a) (P u ,u) plane. If c 2 k + c u > (case 1), the 
motion in the (P u ,u) plane is libration. If c 2 k + c u < (case 2), the motion is rotation, b) (P v ,v) 
plane. Vertical lines v = icosh -1 a/c correspond to a potential wall of the billiard. The region 
between the lines corresponds to phase points moving in the billiard. Shaded region outside the 
lines corresponds to phase points outside the billiard (we do not consider their motion). In the case 
1 (c?k + c v > 0) a phase point which begin to move from a point on the "left" wall with P v > goes 
to the right until a collision with the "right" wall, then jumps down instantaneously to the point 
located symmetrically about the line P v = 0, then moves to the left until collision with the "left" 
wall, then jumps to the initial point, and so on. In the case 2 (c 2 k + c v < 0) a phase point which 
begin to move from a point on the "left" wall with P v > goes downwards until a collision with 
the same wall, then jumps to the point located symmetrically about the line v = 0, then moves 
upwards until a collision with the "right" wall, then jumps to the initial point, and so on. 
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FIG. 2: Jumps of the adiabatic invariant I v on (2,1) resonance. Parameters of the system: e = 
8*10 -5 , uj = 0, di change harmonically with time: d a = d\(l+A\ cos(et)), c4 = (^(l+^b cos(et+0)), 
where A x = 0.3, A 2 = 0.15, d 1 = 2, d 2 = 1, (p = vr/3. 
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FIG. 3: A single jump of the adiabatic invariant on a (2,1) resonance. Parameters are the same 
as in fig 2. 
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FIG. 4: Jumps of adiabatic invariant I v and captures of a phase point into a (2,1) resonance. The 
captured point moves along a resonant curve until it escapes from the resonance. Parameters of 
the system: e = 10~ 3 , u = 2* 10 -4 , d{ change harmonically with time, as in figures 2-3 except that 
A 2 = 0.2 
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FIG. 5: Dynamics of the adiabatic invariant I v of a phase point captured into the resonance (2,1). 
The phase point remains captured forever. Parameters are the same as in Fig. 4 (initial conditions 
are different). In Fig. a) and b) scales are different. 
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FIG. 6: Dynamics of the adiabatic invariant J. Parameters and initial conditions are the same as 
in Fig 5. 
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FIG. 7: Regular dynamics of the system with a/c =const. Coordinates (x,y) of a particle at 
the moments of collisions with a wall are shown. Parameters of the system: e = 0.001, u = 0, 
A\ = 0.15, A2 = 0.15, d\ = 3, c/2 = 1, (f> = 0. a) Time of the integration is one-half of the slow 
period. b)Time of the integration is equal to 2 slow periods. 
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FIG. 8: Irregular dynamics of the system with a/c /const. Coordinates (x,y) of a particle at 
the moments of collisions with a wall are shown. Parameters of the system: e = 0.001, u = 0, 
A\ = 0.15, A2 = 0.25, d± = 3, di = 1, (f> = 7r/3. a) Time of the integration is one-half of the slow 
period. b)Time of the integration is equal to 2 slow periods. 
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